Switching of the vortex polarity in a magnetic nanodisk by a DC current 
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We study the dynamics of a vortex state nanodisk due to a dc spin current, perpendicular to 
the disk plane. The irreversible switching of the vortex polarity takes place above some threshold 
current. The detailed description of these processes is obtained by spin-lattice simulations. 
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The spin torque effect , which is the change of magne- 
tization due to the interaction with an electrical current. 



Heisenberg exchange and dipolar interactions 



was predicted by Slonczewski [T] and Berger |5] in 1996 



During the last decade this effect was tested in differ- 
ent magnetic systems [21 |4l|5] and nowadays it plays an 
important role in spintronics [5J |7]. Recently the spin 
torque effect was observed in vortex state nanoparticles. 
In particular, circular vortex motion can be excited by an 
AC |S] or a DC [H] spin-polarized current. Very recently 
it was predicted theoretically [TO] and observed experi- 
mentally [llj that the vortex polarity can be controlled 
using a spin-polarized current. This opens up the possi- 
bility of realizing electrically controlled magnetic devices, 
changing the direction of modern spintronics ^12] . 



It was shown in [TUl that in easy-plane Heisenberg 
magnets a spin-current which flows perpendicular to the 
nanoparticle plane acts as an effective DC magnetic field 
making energetically unfavorable one of vortex polarity 
states. In this Letter we study the magnetic vortex dy- 
namics in nanodots excited by the spin-polarized cur- 
rent, using the pillar structure, described in detail in 
Refs. [inillS]. We show that the dipolar interaction cru- 
cially changes the physical picture of the vortex dynam- 
ics, breaking the axial symmetry of the system, i.e. the 
z-component of the momentum can be not conserved. 
Qualitatively speaking the dipolar interaction makes two 
main effects |14j : (i) there appears an effective uniaxial 
anisotropy of the easy-plane type, which is caused by the 
faces surface magnetostatic charges and (ii) there appears 
a nonhomogeneous effective in-plane anisotropy, which is 
caused by the edge surface charges (surface anisotropy). 
Due to the surface anisotropy the magnetization near the 
disk edge is constrained to be tangential to the bound- 
ary, which prevents its precession near the edge. That is 
why a simple picture of rotational vortex, which perfectly 
works for the Heisenberg magnet flO' should be revised 
for the nanodot with account of the dipolar interaction. 



The magnetic energy of nanodots consists of two parts: 
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Here Sn is a unit vector which determines the spin 
direction at the lattice point n, i = ^A/(/ioAf|) is 
the exchange length {A is the exchange constant, /io is 
the vacuum permeability, Ms is the saturation magne- 
tization), the vector S connects nearest neighbors, and 
Emn = {n — m)/|n — m\ is a unit vector. The lattice 
constant is chosen as a unity length. 

It is known that as a result of competition between 
the exchange interaction and the dipolar one the ground 
state of a thin magnetically soft nanodisk is a vortex 
state: the magnetization lies in the disk plane XY in the 
main part of the disk being parallel to the disk edge, 
forming the magnetic flux-closure pattern characterized 
by the vorticity q = +1. At the disk center the magneti- 
zation distribution forms a vortex core, which is oriented 
either parallel or antiparallel to the z-axis. The former is 
characterized by a polarity p = -1-1 and the latter p = —1. 
When electrical current is injected in the pillar structure, 
perpendicular to the nanodisk plane, it influences locally 
the spin Sn of the lattice through the spin torque [U |2] 
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Here j — Je/Jp is a normalized spin current, Je 
is the electrical current density, Jp = /ioAf||e|c?/fi, 
d is the disk thickness, e is the electron charge, 
A = 4773/2/ [3(1 + r/)3 - 16r/3/2] , S = (1 + 

77)^/ [3(1 + 77)3 - 1677^/^] , and 7/ e (0; 1) denotes the de- 
gree of the spin polarization; cr = ±1 gives two directions 
of spin-current polarization. 

The spin dynamics of the system is described by the 
modified Landau-Lifshitz-Gilbert equation 
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FIG. 1: (Color online) The time evolution of the vortex switching process from simulations: the top row shows the 3D 
distribution of the magnetization z-component, the bottom row corresponds to the in-plane magnetization distribution near 
the vortex core. Isosurfaces Sx ~ (black curve) and Sy = (orange curve) are plotted to determine the vortex position. The 
current j = —0.1. 



Here the overdot indicates the derivative with respect to 
the dimensionless time t = LUot with ujq = 4ttjAIs, a <C 1 
is a damping coefficient, and !K is the Hamiltonian (jlj. 

To study the vortex dynamics, we have performed nu- 
merical simulations of the discrete spin-lattice Eq. ([s]). 
We consider the case of thin nanodots where the magne- 
tization does not depend on the z-coordinate. We have 
integrated numerically the set of Eqs. ^ over the square 
lattice using the fourth-order Runge-Kutta scheme with 
the time step At = 0.01 and free boundary conditions. 
The lattice is bounded by a circle of diameter 2L. In 
most of the simulations 2L = 100, h = 10, £ = 2.65, 
a = 0.01, cr = -Hi and 77 = 0.26. 

As an initial condition we use the vortex centered in 
the disk origin, see Fig. [T]^a). To identify precisely the 
vortex position we use, similar to Ref. [16j . the crossec- 
tion of isosurfaces Sx = and Sy = 0, see the bottom 
row of Fig. [1] The vortex dynamics results from the force 
balance between a driving force (by the current) , a dipo- 
lar force, a gyroscopical force, and a dissipative force. In 
the simulations we observe that when the vortex and the 
spin-current have the same polarization {jap > 0) the 
vortex does not quit the center of the disk. However, for 
jap < (p = -1-1, cr = -1-1 and j < in our case) the vor- 
tex under the action of the current starts to move out of 
the disk center following a spiral trajectory, see Fig.|2ja). 
The spiral type of motion is caused by the gyroscopical 
force, which acts on a moving vortex perpendicular to 
its velocity in the same way as a Lorentz force acts on 
a charged particle in a magnetic field. The role of the 
charge plays a 7T2 topological charge Q — qp/2. The 
sign of Q determines the direction of the vortex motion, 
which is clockwise for p = +1. At some point (marked 
on Figs. [2] by the green symbol) the vortex switches its 
polarity (p — 1). As is seen from Fig. [l] the mecha- 
nism of the vortex switching is very similar to the one 
observed in other systems [lU [T71 [TSl [111 UHl I2I]- The 



moving vortex excites a non-symmetric magnon mode 
with a dip situated towards the disk center. When the 
vortex moves away from the center, the amplitude of the 
dip increases, see Fig. ^h). When the depth of the dip 
reaches a minimum (Sz — —1), a pair of a new vortex 
and antivortex is created, see Fig. [ijc). The reason why 
the new-born vortex tears off his partner has a topologi- 
cal origin. The gyroscopic force depends on the the total 
topological charge Q. Therefore it produces a clockwise 
motion for the original vortex (q = l,p = 1,Q = 1/2) and 
the new-born antivortex {q = —l,p= — 1,Q = 1/2) while 
the new-born vortex (q = l,p = —l,Q = —1/2) moves 
in the anti- clockwise direction. As a result the new vor- 
tex separates from the vortex-antivortex pair and rapidly 
moves to the origin, see Figs, [ije) and|2j^b). The attrac- 
tive force between the original vortex (q — 1) and the 
antivortex (q = — 1) facilitates a binding and subsequent 
annihilation of the vortex-antivortex pair, see Fig.[ljd). 

The switching process has a threshold behavior. It oc- 
curs when the current \j\ > jsw, which is about 0.012 
in our simulations, see Fig. [3j For stronger currents the 
switching time rapidly decreases. Using typical parame- 
ters for permalloy disks [10l[2T] (r? = 0.26, A = 26 pJ/m, 
Ms = 860 kA/m, a = 0.01), we estimate that the time 
unit 1 /loq = 50 ps, the critical current density is about 
0.1 A//im^ for a nanodot of 20nm thickness. The total 
current for a disk with diameter 200 nm is about 10 niA. 

To summarize, we have studied the magnetic vortex 
switching under the action of a DC electrical current. 
We showed that the switching mechanism is essentially 
the same as in the cases when it is induced by a magnetic 
field pulse [HI El [H] , by an AC oscillating ^ or ro- 
tating field [21], or by an in-plane electrical current [TT| . 
There are two key points in this process: (i) the dipolar 
interaction causes the deformation of the magnetization 
profile for the fast moving vortex, which finally results 
in the creation of an additional vortex-antivortex pair. 
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FIG. 2: (Color online). The vortex dynamics for j — —0.1. The vortex trajectory before the switching (a) and after it (b). The 
radius of the vortex trajectory as a function of time (c): at r = 877 the vortex polarity is switched. 
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FIG. 3: (Color online). Switching time as a function of the 
applied current. 



(ii) the topological charge structure of these three exci- 
tations secures the survival of the vortex with the new 
polarity direction or, in other words, the polarity switch- 
ing. The detailed study of the vortex dynamics including 
the switching process is under construction. 
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